Clock shifts of optical transitions in ultracold atomic gases 
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We calculate the shift, due to interatomic interactions, of an optical transition in an atomic Fermi 
gas trapped in an optical lattice, as in recent experiments of Campbell et ai, Science 324, 360 
(2009). Using a pseudospin formalism to describe the density matrix of the internal two states of 
the optical transition, we derive a Bloch equation which incorporates both the spatial inhomogeneity 
of the probe laser field and the interatomic interactions. Expressions are given for the frequency 
shift as a function of the pulse duration, detuning of the probe laser, and the spatial dependence of 
the electric field of the probe beam. In the low temperature semiclassical regime, we find that the 
magnitude of the shift is proportional to the temperature. 
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In the continuing quest to develop improved atomic 
clocks, attention has turned to fermionic atoms because 
shifts of transition frequencies due to interatomic interac- 
tions (so-called clock shifts) are expected to be strongly 
suppressed in gases of identical fermions due to the Pauli 
exclusion principle. It came as a surprise that nonzero 
interaction shifts were observed in a recent experiment 
on the 1 So~ 3 Po optical transition for 87 Sr atoms trapped 
in an optical lattice P, 0], since for a gas of fermions 
it has been shown theoretically that for a homogeneous 
probe field there should be no frequency shift [3j. Thus 
inhomogeneity is indispensable for observing nonzero fre- 
quency shifts in identical fermion samples. The authors 
of Ref. [l| attributed the shifts to the combined influence 
of the spatial inhomogeneity of the probe laser field and 
the interatomic interaction. The same underlying phys- 
ical mechanism gives rise to the Leggett-Rice effect in 
spin transport in liquid 3 He jj], and spin segregation in 
ultracold bosons Q and fermions 



In this paper, we introduce pseudospin operators to 
describe atomic correlations within each motional states 
of the atoms and derive a Bloch equation that describes 
the evolution of the pseudospin under the combined effect 
of an external probe laser field and interatomic interac- 
tions and derive expressions for the frequency shifts to be 
expected under the experimental conditions of Ref. 
Since this work was largely completed, we became aware 
of Refs. 0, Q in which, for a particular form of the inho- 
mogeneity of the probe field, the problem is approached 
using the wave function of the state. 

Basic formalism. In the experiment, atoms of 87 Sr, a 
fermionic isotope, were initially prepared in one hypcrfinc 
state of the 3 Po excited-state manifold (denoted by e) and 
transferred to a hyperfine state of the 1 So ground-state 
manifold (denoted by g) by application of a probe laser 
field [H- In the absence of the laser field, the system is 



described by the Hamiltonian 



H = I d 3 r 
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where the ip a are annihilation operators for atoms in 
the internal states |a) and E a are the internal energies. 
The effective low energy interaction coupling is given by 
g = Ai:a s /m (K = 1) where a s is the s-wave scattering 
length for collisions between a ground state atom and 
an excited state one. The external dipole optical poten- 
tial V(r) is a superposition of a deep one-dimensional 
optical lattice potential aligned along the z direction for 
which the frequency of oscillations about the minima of 
the potential is ^lattice ~ 80 kHz, and a harmonic po- 
tential with uj± — uj x — Lu y w 450 Hz and u z <C 80 kHz. 
(The potential V(r) is the same for the two atomic states 
since the wavelength of the optical lattice is chosen to be 
Aiattice ~ 813.429 nm, the magic wavelength at which 
the electric polarizabilities are the same for the ground 
and the excited states.) Atoms are thereby confined to 
a number of pancake-shaped regions extended in the x- 
and y-directions and centered on the z-axis, and are es- 
sentially limited to the ground state for motion in the 
z-direction under experimental conditions. 

It is convenient to introduce the pseudospin density op- 
erators S(r) = [tpl(r),vj^(r)]cr[vj e (r),ipg(r)] T /2, where a 
are the Pauli matrices. The coupling between the probe 
laser field E(r, t) and the atoms is due to dipole transi- 
tions between the two states, and can be described by 
the Hamiltonian 



^probe 



d 3 r((e|d-E(r,t)| 5 ) ^(r)^(r) + h.c. 



- d- i r[B x (r,t)S x (r)+B y (r,t)S y (r)) (2) 



where B x (r,t) = (e|d • E(r, t)\g) + (g\d ■ E(r,t)|e) and 
B y (r,t) = i({e\d • E(r,t)\g) - (g|d • E(r,*)|e» are the 
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components of the "pseudomagnetic field" Jj|._and d is 
the electric dipole moment. Note that Refs. [7J, [8[ assume 
only B x (r, t) nonzero. However, both B x (r, t) and B y (r, t) 
depend on the details of the probe laser field and are 
generally not zero. 

Block equation. We study the coherent evolution of 
the system governed by H + i? pro bo through the Bloch 
equations for the pseudospin operators. We expand the 
field operators ip a (r) — J2i a m^ , i( r )j where (f>i(r) are the 
free particle eigenstates for H with g = 0, and CLcxi are 
the corresponding annihilation operators. In the frame 
co-rotating with the probe laser, we derive 
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— c.c.) and fli Z = A, where fli 



-J d 3 r|^(r)| 2 (e|d-E_(r)| ff ) and A = E e - E g - lj l is 
the detuning. Equation (|3|) is valid in the weak interac- 
tion limit, g — ► 0, and the Lamb-Dicke regime where the 
recoil of the atoms due to the scattering with the probe 
laser is negligible and it shows that the magnitude of the 
pseudospin Si remains unaltered. 

At this stage, the joint effect of the spatial inhomogene- 
ity of E(r, t) and the interaction is manifest in Eq. ©: if 
all pseudospins Si initially point towards the north pole 
of the Bloch sphere (as in Ref. corresponding to all 
atoms being in the excited state, the first term causes 
the pseudospins of states with different tti to precess by 
different amounts, and then the interatomic interaction 
term no longer vanishes. 

We convert Eq. ([3|) into an integral equation, 

Si(t) =G < (i)S i (0) 

+ 2YV-/ dt' Gi(t — t')Si(t') x Sj(t'), (4) 
j Jo 

where S;(0) are the initial values of the pseudospins. In 
the Cartesian basis (x, y, z), the Green's function satisfies 
the equation 



with 



Q = Q 



dt 
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where Q x 
and f2 = 



f2 sin 8 cos d>, f2„ = Q sin 8 sin 6, ft. 



cos 8 



Q x + fly + Ql . The matrix Q can be diag- 

onalized by the matrix A = (w 1 ,w 2 ,t | 3) where the vec- 
tors are given by vf = (sin 8 cos </>, sin 8 sin (j>, cos 8), v\ = 
(cos 8 cos 4> + i sin 0, cos ^ sin <j> — i cos (j>, — sin 0) / \/2 and 
f -f = (cos cos (/) — i sin cj>, cos sin0 + i cos <j>, — shi8)/^/2. 
Thus 



G(t) = 8(t)A 
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From Eq. ([4|, the deviation of the total pseu- 
dospin S from its value without interaction, S (t) = 
EiGi(*)S?(0),is given by 



where S, = ([a^, a gl \a[a ell a gl ] T )/2 and = 5S (0 

g J d 3 r|0j(r)0j(r)| 2 . Thus the expectation value of the 
total pseudospin is given by S = ^~J j Sj with (...) denot- 
ing the trace with the density matrix. The probe laser 
field has a generic form E(r, t) = E_ (r) exp(— iuijj) + 
El(r) exp(iu>Lt) with ujl the laser frequency. The driv- 
ing fields f2i = (yti x , fliy, Cliz) have components f2j X = 



= VV- / dt' [Gi(t - O - G,(i - t')} [S^t') x S,(t')] , 

(8) 

where the subscript i of G indicates Q. = $7^. Equation 
© shows that the change in the pseudospin is zero if S~2^ 
is independent of i. 

Frequency shift. In the experiments the frequencies 
Sttij = (S~2 j — flj)/2 and g^ are small compared with 
the average tti and therefore we may treat them per- 
turbatively. For the interatomic interaction this means 
that the change in the total pseudospin is a sum over 
contributions from independent pairs of motional states, 
and therefore we may consider just two motional states 
(which we shall denote by 1 and 2), and then sum over 
all possible pairs at the end. From Eq. ([3|), the total 
pseudospin of the pair of states S+ = Si + S2 and the 
difference S_ = Si — S2 satisfy the equations 



dt 



dt 



S+ =fi x S+ +5ft x S_ 
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S_ =n x S_ +5ft x S+ - 2g 12 S + x S_, (10) 



with ft = (Six + n 2 )/2 and Sfl = (fix - fl 2 )/2. Note 
the interaction does not appear in the equation for S + 
explicitly. This is consistent with the argument based on 
the SU(2) symmetry of the interaction [3j]: ii Sfl = 0, the 
dynamics of S + is only governed by l~i and no interaction 
effects can show up. The physical mechanism for creat- 
ing a frequency shift due to interatomic interactions is as 
follows. Initially all pseudospins are aligned in the same 
direction, and therefore S+ x S_ is zero and the inter- 
action has no effect. However, according to Eq. (fTUf the 
field inhomogeneity creates a component of S_ perpen- 
dicular to S + , and therefore the interactions can have an 
effect. Subsequently, according to Eq. |9]) this results in 
a change in S+. Consequently, the leading contribution 
to the shift is proportional to g(8fl) 2 . 
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For the initial conditions of interest, Si(0) = Si(0)e 2 
and S 2 (0) = S 2 (0)e 2 , where $(0) = fr/2 and is the 
initial distribution function. Solving for the change of 
S+(f) due to the interaction to second order in dfl and 
first order in g\i by iterating integral equations similar 
to Eq. (0]), we obtain 

SS+(t) 

=2 / dt' [ dt" [ dt"'G(t - t')SSl x G(t' - t") 
Jo Jo Jo 

ffia{[G(t")S-(0)] x [G(i"-i"')mxG(t"')S-(0)] 

+ [G(t" - t"')Sn x G(i"')S+(0)] x [G(t")S+(0)] }. 

(11) 
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The processes corresponding to this equation are that at 
time t'" the diffence in precession frequencies gives rise to 
a difference in orientations of the two pseudospins, then 
the interatomic interaction acts at time t" and finally the 
difference of precession frequencies leads to a change in 
the total pseudospin of the pair of states. After simplifi- 
cation Eq. (fTTj) becomes 



5S+(t) 

t t' t" 

---&9nSx(Q)S 2 (0) I dt' [ dt" f dt'"G(t-t')Sfl 
Jo Jo Jo 

x G(t' - t")\G(t")e z x [G{t" - t'")5U x G(t"')e z ]\. 

(12) 



The quantity measured in experiment is the population 
of atoms in the ground state, P g — N/2 — £\ S{ Z , where 
N is the total number of atoms. On summing Eq. (TT2|) 
over all pairs of motional states one finds 



6Pg(t) = ClHl + C 3 H 3 + C 3 S 3 



(13) 
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FIG. 1: Contour plot o f the frequency shift Sui, in units of 
(2N/Q%.y) SijiP'i ~ ^j) 2 as a function of dimensionless detun- 
ing, 5 = A/Q. xy and dimensionless pulse duration, r = Q xy t p . 
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2(l + 5 2 ) 7/2 
x cos Fr\/l + S 2 ^ - \J\ + S 2 sin 2 (iVl + 6 2 ) 

-r (-1 + S 2 + 25 A ) sin (iV 1 + S 2 )} , (18) 
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Ae z , SSlij\i 



Here 



with 8 = A/tt xy , t = Q xy t, Vt xy = fl t - 
(SSlij ■ fl X y)fl xy /fll y , and Sflijx = Sfl t 
the bar denotes an average over motional states. 

Experimentally, the pulse duration t p is fixed and pairs 
of detunings Ai and A 2 which render the same popu- 
lation P g in the ground state at the end of the pulse, 
i.e., S z (t p , Ai) — S z (t p ,A 2 ), are measured. The clock 
shift is defined as 8u> = (Ai + A 2 )/2. The reason why 
this method can detect the effects of interatomic inter- 
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actions lies in the symmetry of the Bloch equation ([3]). gij ~ [min(n,- x , nj X ) mm(m y , nj y )] x l 2 . The difference in 
For simplicity, let us assume f2j a = 0. Without inter- the driving frequencies (fij — J~ij) 2 ~ in%x + — — 



actions, Eq. ([3]) is invariant under the transformation 
A, -S^S^S^}, which implies 
Ai = — A 2 and &j = 0. When g ^ 0, Eq. |S]) becomes in- 
variant under the transformation {A, g, Si X , Si y , S{ Z } — > 
{— A, — — Ste, Si B , 5i Z }, which indicates &t> 7^ if fij 
are different. To order g(Stt) 2 , from Eq. (fl3|) . we have 

2^ z (t p ,A) 



5w 



(20) 



dS°(t p ,A)/dA' 

Field inhomogeneity. To compare our results with ex- 
periment we need to take into account the form of the 
probe laser field, which is a linearly polarized Gaussian 
beam. With the assumption that the axis of symmetry 
of the laser field coincides with that of the optical trap, 
the electric field vector is given by 



E_(r) 



- «>(0) 
=E Q e x ——exp 



x 



y 



x exp 



w(z) ^ \ w 2 {z) 
2tt Tr(x 2 +y 2 ) 
X L X L R(z) 



»C(*) , (21) 



with the field width w(z) = w(0)y/l + (z/zr) 2 , the ra- 
dius of curvature R(z) = z [l + {zr/ 'z) 2 ~\ and the Gouy 
phase C(z) = arctan(z/z^). The Rayleigh range is 



zr — t:w(0) / \l- The laser wavelength is Xl 
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nm and the beam divergence is estimated to be 9 = 
X l /ttw{0) ps 0.01 Q. Atoms interact with each other 
only within a single pancake. For motional states spec- 
ified by the quantum numbers for the three Cartesian 
directions = (n ix , n iy , 0) associated with a given pan- 
cake, the spread of the directions of S7j over the harmonic 
motional states is smaller than 10 -3 in the temperature 
regime T ~ 1 [iK. Therefore we can neglect terms in- 
volving S£l± in Eq. (fT3| . In Fig. JT]) we show contours 
of constant frequency shift as a function of pulse dura- 
tion and laser detuning for conditions of experimental 
interest. In agreement with experiment, we find that the 
shift can vanish, even thougho is nonzero. However, con- 
trary to initial expectations [l| , vanishing of the shift does 
not correspond to there being equal numbers of atoms in 
the ground and excited states at the end of the pulse, 
P g (t p ) = N/2. From Fig. (JTJ), for r < ir, we see that for 
longer pulses, the zero shift occurs at smaller detunings, 
i.e., larger P ff . As the temperature decreases, atoms will 
concentrate more around the center of the trap, leading 
to a larger Q xy , e.g., longer pulses fl X yt P (for fixed t p ). 
Thus we predict that as the temperature is lowered, the 
zero shift point will move to larger P g , a result consistent 
with experiment [l(. 

The temperature dependence of the magnitude of the 
frequency shift can be easily extracted in the regime 
< T < 1/tow 2 (0). The coarse-grained den- 
sity distributions associated with the two-dimensional 

harmonic oscillator wave functions have the semiclassi- 

r /2 

cal form ~ [(n x — mujx 2 /2){n y — mwx 2 /2)] . Thus 



nj y ) . Since typical values of the n x and n y are propor- 
tional to T, the average of the prefactor Si in Eq. (fl3|) 
<7y(fij — ^j) 2 behaves as T. The ratio between the mag- 
nitude of the shift measured at 3 uK and that at 1 /xK is 
about three [lj], in agreement with this behavior. 

In this paper we have derived a closed expression for 
the dependence of the clock shift on the duration of the 
probe pulse, the laser detuning, and the geometry of the 
probe laser beam. In the calculations we have taken into 
account the fact that the effective exchange interaction 
between atoms in different motional states depends on 
the specific states in question, and that the "Rabi fre- 
quency" fli depends on the motional state. As a check 
on our understanding of the basic physics it would be 
valuable to confirm the predicted dependences. In par- 
ticular, one could make experiments for the case when 
the probe beam is not collinear with the axis of the trap. 

We are very grateful for many valuable discussions with 
Jan Thomsen on the experiments of Ref. [l[ . Useful con- 
versations with Jun Ye are also acknowledged. 



[6] 



G. K. Campbell, M. M. Boyd, J. W. Thomsen, M. J. Mar- 
tin, S. Blatt, M. D. Swallows, T. L. Nicholson, T. Fortier, 

C. W. Oates, S. A. Diddams, N. D. Lemke, P. Naidon, 
P. Julienne, Jun Ye, and A. D. Ludlow, Science 324, 360 
(2009). 

S. Blatt, J. W. Thomsen, G. K. Campbell, A. D. Ludlow, 
M. D. Swallows, M. J. Martin, M. M. Boyd, and Jun Ye, 
larXiv:0906.T419l 

A. J. Leggett, Phys. Rev. Lett. 29, 1227 (1972); Z. Yu and 

G. Baym, Phys. Rev. A 73, 063601 (2006); G. Baym, C. 
J. Pethick, Z. Yu, and M. W. Zwierlein, Phys. Rev. Lett. 
99, 190407 (2007). 

A. J. Leggett and M. J. Rice, Phys. Rev. Lett. 20, 586 
(1968); A. J. Leggett, J. Phys. C 3, 448 (1970); L. R. Cor- 
ruccini, D. D. Osheroff, D. M. Lee, and R. C. Richardson, 
Phys. Rev. Lett. 27, 650 (1971). 

H. J. Lewandowski, D. M. Harber, D. L. Whitaker, and 
E. A. Cornell, Phys. Rev. Lett. 88, 070403 (2002); M. 6. 
Oktel and L. S. Levitov, ibid., 230403 (2002); J. N. Fuchs, 

D. M. Gangardt, and F. Laloe, ibid., 230404 (2002); J. 

E. Williams, T. Nikuni, and C. W. Clark, ibid., 230405 
(2002); 

X. Du, L. Luo, B. Clancy, and J. E. Thomas, Phys. Rev. 
Lett. 101, 150401 (2008); X. Du, Y. Zhang, J. Petricka, 
and J. E. Thomas, Phys. Re v. Lett. 103, 010401 (2009). 
K. Gibble, larXiv:0908.3147l 

A. M. Rey, A. V. Gorshkov, and C. Rubbo, arXiv: 
0907.2245. 

To avoid making the notation too clumsy we denote direc- 
tions in coordinate space and in pseudospin space by x, y, 
and z, but the reader should bear in mind that the two 
spaces are distinct. 



